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I. INTRODUCTION 



The relativistic heavy ion collider (RHIC) at BNL is a unique facility which collides two 
heavy ions at * = 200 GeV to study the production of a —on plasma Q aad 
two polarized protons to study the parton spin structure of the proton [2j. Measurements of 
heavy probes such as J / ip production and Drell-Yan production are useful tools to study the 
quark-gluon plasma in heavy ion collisions and to extract the polarized gluon distribution 
functions inside the proton in polarized p-p collisions Q, 0]. Understanding the correct 
production mechanisms for these processes is important. Heavy quarkonium production can 
be described via the non- relativistic QCD (NRQCD) color-octet mechanism 0,0]. 

In NRQCD the energy eigenstates of heavy quarkonium bound states \H > are labelled 
by the quantum numbers J PC , with an additional superscript to give the color; (1) for 
singlet and (8) for octet. When the Fock states are analysed then the dominant component 
in S-wave orthoquarkonium is the pure quark- ant iquark state |QQ[ 3 >S'j 1 ' ) ] >. A state with 
dynamical gluons, such as \QQ[ 3 Pj 8 ^]g > does contribute but only with a probability of order 
v 2 , where v is the typical velocity of the non- relativistic heavy quark (and antiquark). The 
other states, such as \QQ[ 3 S[ 1,8) ]gg >, IQQ^S^g > and \QQ^D ( } fi) }gg > contribute to the 
probability in even higher orders in v (see later). Correspondingly the dominant states in 
P-wave orthoquarkonia are the states IQQP-Fj 1 ^] >, and the states with dynamical gluons 
such as |<5Q[ 3 5'j 8 ' ) ](7 > contribute with a probability of order v 2 . After a QQ is formed in a 
color octet state it may emit a soft gluon to transform into the singlet state |Q<5[ 3 -Pj 1 ' ) ] > 
and then become a J/ip state by photon decay. The QQ pair in a color octet state can 
also emit two long wavelength gluons and then become a J/ip state. All these low energy 
interactions are negligible and the non-perturbative matrix elements, labelled by the above 
quantum numbers, can be fitted from experiments or can be determined from lattice field 
theory calculations. 

Using the NRQCD color octet mechanism heavy quarkonia production rates have been 
ea.cu.ated for the p-p Tevatron eollider QQ, for the e-p HERA coHider Q, for the e + -e~ 
LEP collider [10] and for fixed target experiments |rjj. Also the recent PHENIX data for 
J/ ip production in unpolarized p-p collisions can be explained by this same mechanism [12] . 

The RHIC offers a wide variety of measurements with respect to J/tp production. They 
involve J/ip production (with and without definite polarizations), in unpolarized p-p, d-Au, 
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Cu-Cu and Au-Au collisions and in polarized p-p collisions. Since the maximum transverse 
momentum of heavy quarkonium that can be measured at the RHIC is around 10 GeV/c, 
the parton fragmentation contribution to heavy quarkonium production will be very small 
and we will neglect it in our study. The main contributions to heavy quarkonium production 



at the RHIC are the parton fusion processes 
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1$. 



The inclusive heavy quarkonium production cross section (summed over quarkonium po- 
larization states) in unpolarized and polarized partonic collisions were calculated in j^, ^| 
and Q respectively. Heavy quarkonium production cross sections with definite po- 
larizations in unpolarized partonic collisions were calculated in Q|. In this paper we 
will study the inclusive rapidity distributions for heavy quarkonium production with specific 
polarizations in polarized p-p collisions. We will evaluate the partonic level cross sections 
for the processes qq, gg — > J /ip{X){ip'{^)) m polarized p-p collisions where A is the helicity 
(polarization) of the heavy quarkonium state. This study should be regarded as a prelimi- 
nary analysis of the leading order (LO) contributions to charmonium production in specific 
polarization states. The factorization ansatz, which separates the short- distance coefficient 
functions, calculable in perturbative QCD, from the long distance matrix elements, which 
are fitted to data, is only proved if the charmonium state is produced at a large transverse 
momentum [6(. Our LO analysis considers only charmonium production in the forward di- 
rection at a finite rapidity. It will be followed later by a next-to-leading (NLO) calculation 
where additional quark and/or gluon radiation is produced in the final state so that the 
charmonium state does have a finite transverse momentum. The reason we need these re- 
sults is that the PHENIX collaboration at the RHIC will measure J/ip and ip' production 
with definite polarizations in polarized p-p collisions at -Js = 200 GeV and 500 GeV jjj. 



Since polarized heavy quarkonium production at the Tevatron energy scale 



17| is not ex- 



plained by the NRQCD color-octet mechanism it will be useful to compare our results 
for J/ip and ip' polarizations with the future data at the RHIC. The study of polarized heavy 
quarkonium production in polarized p-p collisions at the RHIC is also unique in the sense 
that it probes the spin transfer processes in perturbative QCD (pQCD). Note that decays 
from higher quarkonium states to the J/ip are ignored. We only consider direct production. 

The spin projection method is used to evaluate the inclusive cross section for heavy 
quarkonium production (summed over polarization states) in parton fusion processes |8|. 
However, the cross section for heavy quarkonium production with a specific polarization in 
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the final state can involve additional matrix elements that do not contribute when the polar- 
ization is summed. This includes interference terms between partonic processes that produce 
heavy quark-antiquark pairs with different total angular momenta. Such interference terms 
cancel upon summing over polarizations. These interference terms can be calculated by 
using the helicity decomposition method ^(|. Hence we will use the helicity decomposition 
method to calculate the square of the matrix elements for heavy quarkonium production 
with definite helicity in polarized partonic collisions. Using these results we will compute 
the rapidity distributions of the cross sections and spin asymmetries of heavy quarkonium 
production with definite helicity states in polarized p-p collisions at RHIC at y/s = 200 GeV 
and 500 GeV within the PHENIX detector acceptance ranges. 

The paper is organized as follows. In section II we derive the partonic level cross sections 
for heavy quarkonium production with definite polarizations in polarized q-q and g-g parton 
fusion processes using the helicity decomposition method within the NRQCD color-octet 
mechanism. In section III we present the results for the differential rapidity distributions 
and spin asymmetries for the J/ip(X) and in the PHENIX detector acceptance range 

in polarized p-p collisions at ^/s = 200 GeV and 500 GeV. We then discuss these results 
and give our conclusions. 



II. INCLUSIVE HEAVY QUARKONIUM PRODUCTION WITH DEFINITE 
HELICITIES IN POLARIZED PARTONIC COLLISIONS 

In this section we will use the NRQCD color-octet mechanism and derive the square 
of the matrix element for inclusive heavy quarkonium production with definite helicities 
in polarized partonic fusion processes. We will consider the (polarized) partonic fusion 
processes qq — ► H(X) and gg — > H(X) where A is the helicity of the produced heavy 
quarkonium state H. We will use the helicity decomposition method [16] within the NRQCD 
color-octet mechanism to calculate these processes where both initial and final state particles 
are polarized. 



4 



A. The qq fusion process 



The production of a heavy charmonium state with helicity A, where A = 0, ±1 correspond 
to longitudinal and transverse polarization states respectively, begins with the calculation 
of the production of a heavy quark anti-quark pair. The matrix element for the light quark- 
antiquark (qq) fusion process q(ki) + q(k 2 ) — > Q(pi) + Q(p2) producing a heavy quark- 
antiquark (QQ) pair is given by 

M q ^ QQ = ^v^h^u^u^h^vfa), (1) 

where P^ = V\ + P2 = K + K and pt = + Lfq j and p% = P^/2- L^q j . Here P^ 
is the CM momentum of the pair and q l is their relative momentum in the CM frame. The 
latter vector does not have any time component so i — 1, 2, 3 only. Lj is the boost matrix 
defined in with both Lorentz and three vector indices. In terms of non-relativistic heavy 
quark Pauli spinors (£ and 77) we obtain (up to terms linear in q): 

\M g ^ QQ \ 2 = ^ 1,'V'n' if .(fcij^nf^) ,(fc 2 ) 7 /n(^i) ^^'A , (2) 

where m is the mass of the heavy quark. For massless incoming quarks and antiquarks we 
have 

u(ki)u(h) = - (1 + /ii7 5 ) 

v(k 2 )v(k 2 ) = ^(1-/1275)7^2- (3) 

The polarized partonic matrix element squared involves the helicity combination (+, +) — 
(+, — ) with +, — denoting the helicities hi, h 2 of the incoming partons. Then from eq. (0) 
we find 

A\M q ^ QQ \ 2 = l-^^a'T^'^a^rj [2m 2 n inj - S^h-h)}, (4) 

using (k 2 ■ L)i = —(k\ ■ L)i = m rii. Here ni,nj are the components of unit three-vectors 
ni, n 2 which specify the polarizations of the heavy quarks and heavy antiquarks respectively 
in the charmonium bound state. Their ^-components are usually chosen along the beam 
direction. In a frame where P, k 1; and k 2 are collinear, they are also collinear with the third 
components of ni, n 2 . Taking the leading order term in an expansion in q so that P 2 = 4m 2 
we obtain 

A \M q ^ QQ \ 2 = £ [n t n 3 - t/Vt^ eVT^ . (5) 



Averaging over the initial color (by dividing by 9), we therefore get 

A 2 2 

A \M q ^ QQ \ 2 = [n t n 3 - 8 i3 ] . (6) 

The two-component spinor factors can be identified with various heavy quarkonium bound 
states H{X) with different quantum numbers as follows, (see Appendix B in |l6| |) 

4m 2 rj'^'^V = <xHPh { x) ^X> = \rn < Of (%) > , 
4m 2 ^Tr^T a V = < X j T a ^P HW ^T a X > = ~m < O*CS ) > , 
WTyVe'eV'ry = < X VtyP H(A) VV'x> = ^L^ A m < Of ( 3 50 > , 
4m 2 77VTVfV'T°77 ee < x VT^P^ (a) ^Vr X > = ^[/^m < Of f^) > , 
4m 2 g^V^'eV^ ee < x t ( _^m )(7 ^ Pif(A) ^(.'ijn)^ > 
= 4f4£/ iA cT"m < Of ( 3 P ) > , 

4m 2 q n q m ^a i T a ^^a j T a V ee < ^{-tp^T 1 ^ P H (x) ^(~D n )a j T a X > 
= AUl U jX 5 mn m < Of ( 3 P ) > • (7) 

The helicity index A is a vector index in the spherical basis. The spherical basis states and 
cartesian basis states are related by a unitary transformation matrix U\i which satisfies the 
relation 

i 

£CW=5ac (8) 

i 

where n l is along the z-direction. Using the above relations we finally obtain 

A 2 2 

A|M^ (A) | 2 = - [1 - 6 X0 ] < Of CS 1 ) > . (9) 

The polarized quark-antiquark fusion process cross section is given by 

A^h(a) = -6($- Am 2 ) ^ [1 - 6 X0 ] < Of ( 3 5x) > , (10) 
and therefore vanishes for A = 0. 
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B. The gg fusion process 



The matrix element for the gluon fusion process g(kx) + g(k 2 ) — » Qijpi) + Q(P2) after 
including s, t, and u channel Feynman diagrams is given by 

M gg ^ QQ = - g 2 e^hyjih) i(^5 ab + ±d abc T c ) + '-f^T* F^} , (11) 



where 



s~ = fl(Pl) |2^LzA±^ + 2^M^!]„ fe) (12) 

2pi ■ fci 2pi • fc 2 



and 



r (A-^+m) 7- _ 7- & - $ 2 + m) 7 ^ 
lPlJ[ 2pi.An 2 Pl -k 2 

+ ^2 k 2 , -h - k 2 ) lx ]v(p 2 ) . (13) 

The three gluon vertex is denoted by V^ x (ki, k 2 , k 3 ) = [{k x -k 2 ) x g^ + (k 2 -k 3 yg uX + (k 3 - 
ki) v g x ^}. Using various identities among the spinors and boost matrices from the appendix 
A of |l6j and after performing a considerable amount of algebra we find 

+ (L J l)n [P"L% ~ P'L) + 2<T(L ■ k{) 3 - (h - k 2 yL) - (h - fc 2 )"L£] q n S}^r, 



+ [P*L v n - P»Ifl q n e* J V + ^ [P^L) - F"2#] (14) 



and 



Hence 



-- [2<T(L • h)j - (h - k 2 yL) - (h - k 2 yL»} ?a j r] 
m j j 

+ l[Ltm-LW]<r?oS V . (15) 

fit 



= T^ih - k 2 ) x e^ x P p ^ V + [ (L " 3 l)j {P v U n - pi*L v n - 1<T{L ■ h) n ) 
2m z mr 

+ - [L%Lj + L^Lj] + —,{L- h) n [(ki - hTL" + (h - fa)"^]] g»£V»7 , (16) 



which is symmetric under k\ <-> A; 2 , /i <-> ^ because (&i • L)» = — (A; 2 • L)$, and 

^ = ' (L 2m fc 4 l)w ( fc i - k2)xe^ x P p q n ev + - tfLjf] , (17) 

which is antisymmetric under the interchanges above. For an incoming gluon with a helicity 
Ai the square of gluon polarization vector can be written as 



el(h, X 1 )e:\k 1 , Ai) = t + ki^ + h^ _ ^ w |iM_] , (18) 



with a similar result for the square of the second polarization vector e^(/c 2 , A 2 )e* 6 (fc 2 , A 2 ) 
Using the relation 



e^k^kl = 2m 2 e i i k n k L?L» , (19) 



from appendix A of 16] and choosing longitudinally polarized gluons we find that 



A|.V„„_„rt| 2 = - !L e pqr e p ' q ' r n r n r > 



99- 

w rca6c*a'>r r r r Q*u,'v' , rpab rp*ab t t t t J7i*u'f'l forA 



where 



5«6 = + Irf abc T c and F ab = lf abc T c . (21) 

6 2 2 J v ; 

The cross terms between S and F vanish because 

gab p*ab g Q<*ab pab (22) 

due to their symmetry properties. Using various properties of the Lf matrices from appendix 
A of Q| and performing some lengthy algebra we finally obtain, after averaging over the 
initial color (by dividing by 64), 

A \ M gg^QQ\ 2 = ~ -jp W^'^V + ^l( n • + (" • <l') n j'<lj ~ \{ n ' * Q') n j n f 

IK IK 

-(n x q') 3 (n x ^tjV ^Vr? + — v *T a Z'?T a ri + — [(n • g)n^, + (n • g'K'^ 

-~(n • g)(n • - (n x x q)j']rpo i 'T a Z'(*o j T a ri 

27 

+— (n.q)(n.q')^TreT a v]- (23) 
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After identifying the different bound states as given in eq. (J7J) and then using eq. (JHJ) we 
obtain 



A|M 9 ^ (A) | 2 = O?( 1 S )> + ™<OfCS )> 



27 

+ A(^ao - 1)[< Of ( 3 Po) > +^ < Of ( 3 Po) >] + A < W^i) >] • (24) 



Hence the partonic level inclusive cross sections for heavy quarkonium production with 
definite helicities A in polarized g-g collisions are given by 

A^h(A) = - S(s - 4m 2 ) [< OffSo) > +Y< Of? So) > 

+^(^ao - 1) [< Of ( 3 Po) > + f < Of ( S fl.) >] + ^ < Of^i) >] • (25) 

III. RESULTS AND DISCUSSION 

Using the formulae derived above we compute the LO rapidity distributions and spin 
asymmetries for the heavy charmonium systems J/ip and ip' in longitudinally polarized 
proton-proton collisions at RHIC. This provides interesting information on the polarization 
state of these heavy charmonium states. In terms of the heavy quark relative velocity 
v, the non-perturbative matrix elements for 0$( 3 Pj) and 0$( 3 Si) production scale like v 7 
and for (9 8 ( 1 S'o) production scales like v 6 whereas those for Cs( 1 -Pi), Oi( 1 S'o), production 
scale like t> 10 and Oi( 3 Pj) production scale like v 11 0, Q]. Hence we will not include 
the latter contributions as they are expected to be small. In particular this means that 
contributions from processes in which the QQ pair is produced in color-singlet states are not 
included. Folding eqs. (|10|). (J23j) with parton densities we find the following cross sections 
in longitudinally polarized proton-proton collisions 

x (8 X0 -1) <O j J^' ) ( 3 S 1 )> + H Af g ( Xl ,2m) A/,(— ,2m) 

x [^(l-ko)<0 8 W) ( 3 P )> - <0 8 W) (%)>]], (26) 

where A/(x, Q)(Ag(x, Q)) denote the polarized quark (gluon) distribution functions inside 
the proton at the scale Q. The corresponding production cross sections for unpolarized 
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proton-proton collisions are 



n 3 a 2 ri 



27 sm 3 J4m 2 /s £1 



r/ ' L [/ ? (x 1 ,2m)/ g -(^,2m) 



x (1-5ao) <0 8 J/ ^ 3 
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,4m 2 



( Sl)> + 32^( Xl ' 2m )^(— ' 2m ) 



X [A( 1 -^0)<^8 W) ( 3 n)> + <O 8 W) ( 1 5 )>]]. 

m z 3 

The spin asymmetry A_li / (A) is given by the ratio of the above cross sections 



A LL (X) 



da(X) ' 



(27) 



(28) 



unpolarized pp collisions, cm energy=2()() GeV (\=0) 
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FIG. 1: Differential rapidity distributions for J/ip production at y^s = 200 GeV. The solid, dashed 
and dot-dashed lines, each to be multiplied by a scale factor of 200, are for unpolarized p-p collisions 
with J/ijj helicities A = in scenarios 1, 2 and 3 respectively. 

We begin with a discussion of the rapidity distribution for J /if) production. Color octet 
contributions have been obtained from an analysis of charmonium transverse momenta dif- 
ferential cross section data from the Fermilab Tevatron, see j^J and |s| . In particular central 
values are known for < Og^( 3 Si) > and the combination 

M J ^CS^\ 3 P (8) ) =< O^CSq) > +3.5 < O^'fPo) > /m 2 , (29) 

together with reasonable error ranges from both statistical and theoretical uncertainties. 
Measurements of the polarized cross sections should allow one to determine the individual 
values for these contributions. All we can do at present is assume plausible values to guide 
experimental investigations. 

Therefore we choose three scenarios. 
1. < O^^So) > and < OI^^Pq) > /m 2 are roughly equal, so we set both to 0.0087 

10 



unpolarized pp collisions, cm energy=20() GeV (X=l) 




FIG. 2: Differential rapidity distributions for J/ip production at ^fs = 200 GeV. The solid, dashed 
and dot-dashed lines, each to be multiplied by a scale factor of 200, are for unpolarized p-p collisions 
with J /ip helicities A = 1 in scenarios 1, 2 and 3 respectively. 

GeV 3 . 

2. < O^^Sq) > is much larger than < O^^Pq) > /m 2 so we set the former to 0.039 
GeV 3 and the latter to zero. 

3. < O^^So) > is much smaller than < O J 8 ^( 3 P ) > /m 2 so we set the former to zero 
and the latter to 0.01125 GeV 3 . 

In all three cases M^ 1 ^ 8 } , 3 P (8) ) « 0.039 GeV 3 , which is the central value in |2o| for 



GRV leading order (LO) parton densities |2lj. We also take < O^^Sx) >= 0.0112 GeV 3 
in agreement with the value in [^J. Note that these values can still vary up and down by 
approximately fifty percent. 

The central arm (forward arm) electron (muon) detector at the PHENIX experiment 
covers the J/ijj rapidity range — 0.5 < y < 0.5 (1 < \y\ < 2). We will present our 
differential rapidity distributions and spin asymmetries for J/ip and ip' production with 
helicities A = 1 and in unpolarized and polarized p-p collisions at ^fs = 200 GeV and 500 
GeV in the above detector acceptance ranges. 

We take the charm quark mass m=1.5 GeV and the mass factorization scale equal to 2m. 



22 



23| and 



We choose 



Several groups have produced polarized parton density sets 

the GRV unpolarized LO parton densities j3] and the GRSV [23| polarized densities. The 
latter authors have a standard scenario and a valence scenario. For simplicity we choose the 
former. Therefore we always use the LO four flavour sets (for the u, d, s and g partons) and 
we set 71/ = 4 in the one-loop running coupling constant and the parton densities. For both 
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polarized pp collisions, cm energy = 200 GeV (^.=0) 
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FIG. 3: Differential rapidity distributions for J/ip production at \/s = 200 GeV. The solid, dashed 
and dot-dashed lines are for polarized p-p collisions with J/ip helicities A = in scenarios 1, 2 and 
3 respectively. 

polarized pp collisions, cm energy = 200 GeV (X=l) 
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FIG. 4: Differential rapidity distributions for J/ip production at y/s = 200 GeV. The solid, dashed 
and dot-dashed lines are for polarized p-p collisions with J/ip helicities A = 1 in scenarios 1, 2 and 
3 respectively. 

parton density sets we use A^ = 175 MeV, so that (mz) = 0.121 at the mass of the Z. 

We first choose helicity A = which eliminates contributions from < Og^( 3 S\) > in eqns. 
f!26|) and (}2~Tj) . Note that this contribution only appears in the quark-antiquark channel. In 
Fig. 1 we present the rapidity differential distributions for J/ip production with A = in 
unpolarized p-p collisions at \fs = 200 GeV. The solid, dashed and dot-dashed lines, each to 
be multiplied by 200, are the differential distributions for the scenarios 1, 2 and 3 respectively. 
As expected there is very little difference between them. Note that the y axis is taken from 
-1 to +3 as in the case of the polarized plots, which will be shown shortly. Fig. 2 contains the 
corresponding distributions when A = 1. Now < Oq^( 3 Si) > contributes to eqn. (|2*7j) and 
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unpolarized pp collisions, cm energy=5()() GeV (A.=0) 
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FIG. 5: Differential rapidity distributions for J/ip production at ^fs = 500 GeV. The solid, dashed 
and dot-dashed lines, each to be multiplied by a scale factor of 1000, are for unpolarized p-p 
collisions with J/tp helicities A = in scenarios 1, 2 and 3 respectively. 

is responsible for the tiny difference between the results for the second scenarios (the dashed 
lines). From this we conclude that the gluon-gluon contributions completely dominate the 
quark- ant iquark contribution for this cm energy. The larger differences between the other 
results are due to the changes in the prefactors multiplying the octet contributions in eqn. 

dZZJ). 

unpolarized pp collisions, cm energy=5()() GeV (X=l) 
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FIG. 6: Differential rapidity distributions for J/ip production at y^s = 500 GeV. The solid, dashed 
and dot-dashed lines, each to be multiplied by a scale factor of 1000, are for unpolarized p-p 
collisions with J ftp helicities A = 1 in scenarios 1, 2 and 3 respectively. 

We now give the corresponding results for longitudinally polarized p-p collisions at the 
same cm energy using eqn. ([26)1 . Fig. 3 has A = and Fig. 4 has A = 1. In both cases 
scenario 2, where < Og^( 3 P ) > /m? = 0, yields, as expected from the signs in eqn. (J2BJ), 
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negative results. We conclude that the polarized results are at least a factor of 200 lower 
then the unpolarized ones but they could be much smaller due to cancellations between the 
color octet contributions in eqn. ()26jl and could possibly be negative. 

polarized pp collisions, cm energy=500 GeV (^.=0) 
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y 

FIG. 7: Differential rapidity distributions for J/iJj production at yfs = 500 GeV. The solid, dashed 
and dot-dashed lines are for polarized p-p collisions with J/ip helicities A = in scenarios 1, 2 and 
3 respectively. 

In Figs. 5 - 8 we repeat the plots in Figs. 1-4 but for y/s = 500 GeV. This time we 
divide the unpolarized values by 1000 so that they fit on the same scales as the polarized 
ones. Also we choose the y-axis from -1 to +2 so that the polarized and unpolarized plots 
can be easily compared. The comments we made earlier about Figs. 1-4 are also valid for 
these plots. Here we see that that the polarized values are at least a factor of 1000 times 
smaller than the unpolarized ones, which is due to the rapid rise of the unpolarized gluon 
density at small x, but again they could be even smaller due to cancellations between the 
color octet contributions. 

For completeness we give in Fig. 9 the rapidity distributions for the spin asymmetry All 
defined in eq. (|2*Hj) for J ftp production in p-p collisions at y/s = 200 GeV. This is for the 
scenario 1. The solid and dotted lines are for helicities A = 1 and respectively. We see 
that the ratios are rather flat over the central rapidity range. For scenario 2 All is negative 
and roughly the same in absolute value. For scenario 3 All is positive and about the same. 
Therefore we do not show the latter two cases. 

Finally Fig. 10 contains the rapidity distributions of the spin asymmetries All for J/ip 
production in p-p collisions at y/s = 500 GeV. The solid and dotted lines are the spin 
asymmetries for helicities A = 1 and respectively. The ratio is also rather flat and smaller 
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polarized pp collisions, cm energy=50Q GeV (A^l) 
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FIG. 8: Differential rapidity distributions for J/ijj production at = 500 GeV. The solid, dashed 
and dot-dashed lines are for polarized p-p collisions with J/ip helicities A = 1 in scenarios 1, 2 and 
3 respectively. 

0.005 
0.0045 

0.004 
0.0035 

0.003 
^ 0.0025 

0.002 
0.0015 

0.001 
0.0005 

°-3 -2-10123 

y 

FIG. 9: Differential rapidity distributions for All at = 200 GeV in scenario 1. The solid and 
dotted lines are for helicities A = 1 and respectively. 

than in Fig. 9, which is due to the increase in the unpolarized gluon density at small x. We 
do not give the plots for scenarios 2 and 3 as they are similar in absolute magnitude. 

Let us now consider ip f production. Central values for < 0$ ( 3 Si) > and the combination 

M^'CS®, 3 P (8) ) =< Of'CSo) > +3.5 < Of ( 3 P ) > /m 2 , (30) 

together with reasonable error ranges from both statistical and theoretical uncertainties are 
given in [20j. However we note that only the former, equal to 0.0046 GeV 3 , is derived from 
data analysis while the latter is not. The data were simply not good enough so it was 
assumed that the ratios of the above color combinations are the same for the J/ip and ip' 
respectively. However we have seen that the contribution from < 0$ ( 3 Si) > in the quark- 
antiquark channel is completely negligible compared to the other color octet contributions 
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from the gluon-gluon channel in this energy range. Hence we only need to look at the ratios 
in to determine the reduction factor for the rapidity distributions for ip' production. 
This is 3.5 so all plots given above in Figs. 1-8 can be used for ip' production by simply 
dividing them by this number. The All ratio plots in Figs. 9 and 10 are not affected since 
this factor cancels between the numerator and denominator. 

In this paper we have calculated inclusive production cross sections and spin asymmetries 
for heavy quarkonium states with definite helicities in polarized proton-proton collisions 
using the non-relativistic QCD color-octet mechanism. We have presented the LO results for 
J/ijj and ip' rapidity differential distributions with definite helicities in polarized p-p collisions 
at <Js = 200 GeV and 500 GeV at the RHIC within the PHENIX detector acceptance range. 
One can see from the figures that the contributions from the A = ±1 states dominate. This 
is explained by the fact that the RHIC is a p-p collider so the contribution from the quark- 
antiquark channel is small compared to the contribution from the gluon-gluon channel and 
the coefficient in front of the < Of ( 3 Pq) > term in eqn. (}26|) is larger for A = 1 than for 
A = 0. The PHENIX experiment should be able to measure these spin asymmetries. The 
study of heavy quarkonium production with definite helicities in polarized p-p collisions is 
unique because it tests the spin transfer processes in perturbative QCD. As Tevatron data 
for heavy quarkonium polarization is not explained by the color octet mechanism Q|, it 
will be interesting to compare these theoretical results with the future experimental data at 
RHIC to test the NRQCD color octet heavy quarkonium production mechanism with respect 
to polarization. A calculation of the px distribution in heavy quarkonium production with 
definite polarization states in polarized p-p collisions at RHIC should also show interesting 
features. 
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